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Nickol'skii  (1959)  ars  used.   The  effects  of  ultraviolet  radiation  due  to 
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CHAPTER  I 
INTRODUCTION 

Two  recent  investigations  have  stirred  again  the  controversy  concerning 
the  possible  presence  of  obscuring  matter  in  globular  clusters.   First.,  Idlis 
and  Nikol'skii  (1959)  reported  sighting  tiny  globule-like  objects  in  the  central 
region  of  the  globular  cluster  closest  to  the  sun,,  Messier  kc      Roberts  (i960) 
investigated  obscuring.,  cloud-like  objects  which  -were"  discovered  in  12  out  of 
32  high  altitude  globular  clusters , 

The  arguments  presented  by  these  investigators  that  the  obscuring  material 
is  an  integral  part  of  the  globular  clusters  seem  convincing,   'They  both  showed 
that  the  probability  of  such  vacant  or  obscured  regions  appearing  as  a  result 
of  statistical  fluctuations  in  stellar  distribution  of  the  cluster  is  extremely 
remote,  if  not  impossible.   Roberts  also  demonstrated  that  for  such  obscura- 
tion to  be  caused  by  intervening  halo  clouds  would  require  a  space  density  of 
such  clouds  much  higher  than  in  the  solar  neighborhood  -  a  result  again  extremely 
unlikely.   It  was  also  noted  that  clouds  or  obscurations  do  not  appear  against 
the  background  of  the  elliptical  galaxies. 

In  this  investigation  then^  we  start  with  the  assumption  that  these  objects 
are  indeed  obscuring  medium  in  the  clusters 3    to  be  analyzed  concerning  both 
their  reaction  with  the  galactic  disk  and  the  possible  environment  within  the 
cluster  which  affects  their  formation  and  behavior.   From  the  investigation 
of  Roberts _,  we  sense  a  possibility  that.,  since  not  all  the  clusters  show  obscur- 
ing matter }    there  is  a  selection  effect  operating  to  make  conditions  for  forma- 
t/ion of  such  objects  and  their  retention  and  persistence  within  the  clusters 
more  favorable  in  certain  cases  than  in  others.   This  approach  then  suggests  the 
possibility  that  there  may  be  small  but  significant  differences  between  individual 
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clusters.   However.,  the  globular  clusters  seem  to  be  a  remarkably  homogeneous 

roup  of  objects  with  respect  to  certain  properties.   Von  Hoerner's  (1957)  study 
of  globular  clusters  showed  that  their  standard  deviation  in  luminosity  is  less 
than  1  absolute  magnitude,,  that  their  deviation  of  diameters  is  approximately 
2jo   of  the  mean  diameter.,  and  that  their  internal  structures  are  very  similaro 
This  homogeneity  suggests  that  we  adopt  a  standard  cluster  model  as  the  frame- 
work within  which  to  make  our  order  of  magnitude  calculations  concerning  the 
interstellar  material,,   Because  of  the  great  amount  of  observation  and  study 
devoted  to  the  cluster  Messier  3^  we  have  generally  adopted  this  object  as  our 
model. 

Of  the  many  basic  assumptions  that  must  be  made  in  an  investigation  of 
this  sort,,  perhaps  the  most  difficult  is  that  concerning  the  physical  state 
and  concentration  of  matter  in  the  galactic  disk,,  lying  outside  the  discrete 
clouds  and  galactic  spiral  arms.   Results  from,  the  study  of  optical,,  inter - 
stellar  absorption  lines.,  emission  nebulae,,  absorption  and  obscuration  of 
starlight,  and  from  radio  astronomy  work  in  both  the  continuum  and  with  the 
21 -cm  line  have  given  a  fairly  consistent  picture  of  the  galactic  structure 
in  the  dense  regions.   However,,  a  search  of  the  literature  reveals  that  con- 
fusion and  controversy  presently  exist  on  the  state  and  density  of  matter 
between  the  arms. 

In  his  study  of  absorption  lines  in  the  optical  region,  Stromgren  (19^8) 
implies  that  the  non-cloud  region  is  permeated  with  a  HI  medium  of  approxi- 
mately  0.1  H  atoms  cm  ^0   Spitzer  (19^9)  suggests  at  the  first  symposium  on 
cosmical  gas  dynamics  that  the  interstellar  medium  is  in  pressure  equilibrium 
except  for  the  relatively  small  volume  of  expanding  H  II  regions  near  hot  stars, 
s  model  would  require  the  diffuse  intercloud  medium  to  have  a  kinetic  tempera- 
ture of  10,000  K,  signifying  that  it  is  primarily  ionized  or  HII.   A  summary 
of  his  model  shows: 
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However t   later  Spitzer  (195M  states  that  the  "atoms  in  this  tenuous  intercloud 

medium  could  be  mostly  neutral  and  still  at  a  relatively  high  temperature",  Dufay  (1955) 
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says  that  regions  of  diffuse  hydrogen  of  10  cm   concentration  could  be  HI  regions 

of  temperature  1000  K  or  less,  because  of  the  long  times  (10  years)  required  to 
reach  equilibrium  at  this  density,   A  gas  at  this  density  cools  at  a  very  slow 

However ,  on  page  1^3 he  remarks  that  much  of  the  intercloud  medium  could 
be  HII.   Finally,  the  question  of  whether  the  medium  is  predominantly  HI  or  HII 
is  discussed  at  length  in  the  Third  Symposium  on  Cosmical  Gas  Dynamics  held  in 
1958»   Spitzer  (1958)  concludes  that  "in  any  case,  this  intercloud  medium  is 
bably  at  a  high  temperature,  so  that,  its  pressure  is  equal  to  that  in  the 
ser,  cooler  HI  regions.   Whether  this  rarified  gas  between  clouds ,  and  not 
near  an  early  type  star,,  is  HI  or  HII  is  an  open  question". 

In  all  probability,  the  intercloud  hydrogen  of  the  galactic  disk  is  not  uni- 
form in  concentration,  temperature,  or  degree  of  ionization.   However,  an  upper 

-3 
limit  to  its  concentration  can  be  set  reasonably  as  0.1  cm   or  one  order  of  magni- 

tude  lower  than  the  concentration  of  1  H  atom  cm   assigned  to  the  spiral  arms  by 
radio  21-cm  line  and  optical,  absorption  line  studies.   Concentrations  as  low  as 
0.01  cm   can  also  be  considered  reasonable. 

In  this  investigation,  we  are  concerned  with  the  gross  dynamics  of  the  encoun- 
ter process  between  a  dense  cloud  or  globule  of  the  globular  cluster  and  a  typical 
cloud  or  the  intercloud  medium  of  the  galactic  disk.   To  simplify  our  calculations 
and  to  enable  us  to  obtain  both  order  of  magnitude  answers  and  limiting  parameters, 
we  have  assumed  the  intercloud  medium  to  be  neutral.   In  fact,  by  using  the 
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macroscopic  approach  of  transfer  of  momentum  and  kinetic  energy  between  large 
regions  of  gas  and  ignoring  the  individual  particle  or  microscopic  behavior ^ 
-we  reason  that  the  question  of  a  neutral  or  an  ionized  condition  for  the  inter- 
cloud  medium.,  interacting  with  the  neutral  clouds  of  the  cluster.,  is  not  a 
sensitive  input  to  the  final  estimates.   Bates  and  Griff ing  (1953)  have  calcu- 
lated collisional  ionization  cross  sections  for  the  neutral -neutral  and  neutral- 
ionized  interaction  of  hydrogen,.   They  find,,  that  at  the  lower  energies  of  approx- 
imately the  range  of  our  problem ^  the  ionization  cross  section  of  the  two  reac- 
tions are  approximately  the  same. 

Little  of  theory  or  experiment 5  applying  to  astrophysieal  problems  of 
a   densities.,  has  been  done  for  the  case  of  plasma  interactions,,   Alfven 
(i960)  has  considered  the  case  of  neutral -ionized  collision.   He  analyzes  the 
case  where  a  non-ionized  gas  collides  with  a  plasma  where  the  kinetic  energy 
per  particle  of  relative  motion  between  the  particles  of  the  two  media  exceeds 
.  ionization  energy  per  particle  of  the  neutral  gas.   He  concludes  that  the 
■  tral  gas  will  be  ionized  and  trapped  in  the  magnetic  field  and  the  density 
of  the  plasma  will  increase.   However.,  implicit  in  his  treatment,,  and  also  in 
the  two  papers  by  Kahn  (1957^  1958)  concerning  collision  of  two  ionized  gases,, 
is  the  assumption  that  the  two  gases  have  comparable  densities.  We  will  find 
that  the  case  which  ultimately  interests  us  in  this  investigation  is  the  colli- 
sion of  a  dense ,   neutral  gas  of  relatively  small  volume  with  a  highly  tenuous 
gas  of  much  larger  extent, 

Kahn  finds  when  two  rarified  masses  of  non-ionized  gas  collide  at  high 
enough  velocity  that  both  gases  will  be  ionized  after  a  relatively  shallow 
penetration,,  and  that  the  counter  streaming  electrons  will  be  stopped  in  a  very 
ort  distance.   This  results  in  a  stationary;,  oscillating  electron  gas  through 
which  protons  are  streaming  in  opposite  directions.   This  so-called  collective 
instability  of  the  electrons  and  not  the  individual  collisions  between  particles 
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is  responsible  generally  for  stopping  the  counter -streaming  gases „ 

However,  Kahn  does  assume  equal  density  of  particles  in  each  stream*,   For 
the  case  where  a  wide  difference  exists  in  the  concentrations  of  the  two  counter, 
streams 3    a  momentum  exchange  treatment  is  appropriate  to  determine  the  result- 
ing motion  of  the  arrested,,  now  unified  conglomerate,, 

So  we  will  restrict  our  attention  in  this  problem  to  the  question  of  what 
happens  to  the  cloud  of  a  globular  cluster  with  respect  to  the  cluster's 
standard  of  rest  when  it  impacts  with  the  medium  of  the  galactic  disk*   We  will 
neglect  considerations  of  individual  particle  reaction. 


-  5 


CHAPTER  II 
DISCUSSION  OF  OBSERVATIONAL  DATA  AND  ASSUMPTIONS 

Idlis  and  Nikol'skii  (1959)  detected  small  obscurations  in  the  central 
regions  of  Uk.      These  globules  appear  twice  as  numerous  photographed  in  the 
blue  as  they  do  in  the  red,,  thus  following  a  rough  1/X   absorption  law  and  indi- 
cating the  presence  of  dust  particles  characteristic  of  the  plane  of  the  Galaxy* 
They  estimated  the  following  average  properties  for  the  globules;   radius  =  0.02  pc., 
optical  thickness  (t)  =  3?   mass  of  dust  =  0,001  solar  mass  (G) 3    total  mass  =  0.  lO,, 
concentration  of  hydrogen  atoms  =  10-cm  ,  total  number  of  clouds  =  1^0.   They 
report  these  globules  as  concentrated  toward  the  center  of  the  cluster  with  maxi- 
mum concentration  at  about  0. 6  pc  distance  from  the  center. 

On  the  other  hand.,  Roberts  (i960)  has  reported  the  detection  of  larger, 
more  massive  objects  in  12  out  of  32  high  latitude  globular  clusters  examined. 
A  typical  cloud  is  that  reported  in  M  13  as  having  a  projected  image  of  1  by  3  pc 
and  a  total  photographic  absorption  of  5  magnitudes.   Roberts  estimates  the  mass 
of  dust  in  this  particular  cloud  as  2, 6  O.   Assuming  the  estimated  itfo   ratio  of 
dust  to  gas  in  the  Galaxy  plane  to  hold  for  the  globular  clusters 3   he  finds  the 
total  mass  to  be  260  O.   Actually.,  for  simplicity  of  calculations.,  this  is  as- 
sumed to  be  the  weight  of  the  atomic  hydrogen  atoms.   This  would  amount  to  a 
hydrogen  atom  concentration  of  the  order  of  lO^cm  .   He  estimates  that  a  typical 
cluster  which  shows  these  clouds  would  have  about  1000  O  tied  up  in  such  objects. 
Roberts  does  not  discuss  the  distribution  or  location  of  his  clouds ,    other  than 
to  show  by  photograph  that  they  generally  appear  at  the  cluster  edge  where  the 
apparent  density  or  luminosity  appear  to  drop  sharply.   However,  from  his  discus - 

)n  of  the  dimensions  of  a  cloud  shown  in  a  photograph  of  M  13,  one  can  deduce 
that  the  cloud  lies  a  minimum  distance  of  about  8  to  10  pc  from  the  cluster  center. 
The  procedure  used  to  determine  the  5  magnitude  absorption  is  not  discussed  by 
Roberts  and  is  probably  a  very  rough  estimate. 
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It  is  clear  that  the  foregoing  two  investigators  are  concerned  each  with 
an  entirely  different  object,  differing  by  ah out  two  orders  of  magnitude  in 
"both  density  and  radius.   Accordingly,  any  discussion  of  the  possible  origins, 
lifetimes,  or  evolutions  of  interstellar  material  in  globular  clusters  should 
concern  itself  and  attempt  to  account  for  both  classes  of  objects.   Throughout 
this  paper,  those  objects  reported  by  Idlis  and  Nikol'skii  will  be  called 
globules  while  those  of  Roberts  will  be  termed  clouds. 

In  making  our  analysis,  we  will  adopt  a  detailed  model  of  a  typical  globu- 
lar cluster.   The  cluster.  Messier  5 s   has  been  observed  and  studied  in  great- 
detail  by  Sandage  ( 195^j  195?)  and  others.   Accordingly,  we  use  this  observa- 

L  data  as  a  basis  for  oar  study.  Also,  we  adopt  Oort  and  van  Herk's  (1959) 
dynamical  model  (d)  of  M3  since  it  corresponds  so  closely  to  its  observed  structure. 

Interstellar  material  could  come  from  any  or  all  of  three  sources:   (l)  resi- 
dual material  left  over  after  original  condensation  of  the  cluster,  (2)  material 
..   has  been  ejected  from  evolving  cluster  stars,  (3)  material  that  has  been 
1 ed  up  from  the  halo  and  galactic  disk  by  the  accretion  of  clouds  within 
the  cluster.   As  we  see  later,  sources  (l)  and  (3)  contribute  a  negligible 
amount  compared  to  (2),   We  adopt  the  results  of  Sandage  for  the  evolution  of 
M3  which  show  that  approximately  100,000  0  have  been  converted  into  interstellar 
material  over  about  10   years  lifetime.   We  also  assume  the  material  to  have 
the  cosmic  proportion  of  elements,  as  adopted  for  interstellar  material  in  the 
galactic  disk. 

Orbits  of  Globular  Clusters-  Most  of  the  rough  orbit  calculations  of  globular 
clusters  are  based  on  the  radial  velocity  measurements  made  by  Mayall  (19^-6)  or 
^oy   Kinman  (1959)°   In  this  investigation  we  will  use  the  more  recent  data  furnished 

f   Kinman,   He  finds  that  the  clusters  follow  highly  eccentric  orbits,  the  average 
eccentricity  being  0,8,   It  is  interesting  to  note  that  a  graph  in  a  recent  paper 
by  Eggen,  Lynden-Bell,  and  Sandage  (1962)  also  shows  an  average  eccentricity  of 

about  0,8  for  high  velocity  stars  of  the  subdwarf  variety,  which  from  ultraviolet 
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excess  measurements,  appear  to  correspond  to  the  main  sequence  of  globular 
clusters. 

Matsunami  and  collaborators  (1959)  have  computed  galactocentric  distances 
for  9^  of  the  118  known  clusters,,  using  Lohmann's  heliocentric  distances*  We  find 
"he  average  galactocentric  distance  to  be  8„5  kiloparsecs  (kpc)  for  their  9k 
clusters.   Since  we  know  that  for  highly  eccentric  orbits^,  an  object  is  most 
apt  to  be  found  near  apogalacticon,  we  make  a  rough  estimate  that  the  average 
apogalacticon  distance  (Q)  is  10  kpc»   Using  the  formula  Q  =  a(l+e),  we  find 
the  semi -major  axis  to  be  about  5°  5  kpc«   Using  the  formula  for  the  latus  rectum 
p  =  a(l-e  )_,  we  find  p  =  2  kpc 

Having  found  the  parameters  for  the  average  orbit,,  we  use  them  to  obtain 

the  average  velocity  (relative  to  the  galactic  frame  of  refernece)  of  a  cluster 

as  it  traverses  the  galactic  disk,,   For  this  calculation  we  assume  that  the 

average  orbital  axis  makes  an  angle  of  V?  degrees  with  the  galactic  disk.,  and 

we  use  the  polar  equation  for  an  ellipse^  r  =  _,  n  ff  ,.  n.  (where  <*  =  anomalv) , 

1+e  cos  <*   x  J ' ' 

solving  for  r  at  the  two  points  of  intersection,,   Upon  finding  r  =  4,.  5  kpc  and 
1.3  kpc ,,  we  select  r  =  3  kpc  as  a  typical  intersection  point  on  the  galactic 
disk. 

We  use  the  Schmidt  model  galaxy  (1956)  for  the  distribution  of  mass  in  the 

P    P 
Galaxy,   We  find  the  following  values  from  his  graph  for  potential  ( in  km  /sec  ) 

in  the  (to,,  z)  plane-  at  Q  =  10  kpc.,  $^  =  30,000  km  /sec  5  at  r  =   3  kpc, 

p    P 
4y  =  d0, 000  km  /sec  „   Adopting  the  value  Vq  =  60  km/sec  at  10  kpc  and  using 

the  formula  ^a*  -   Su  =  Jvr  -  $r,  we  find  V  =  330  km/sec. 

We  have  made  the  above  approximations  for  a  "mean"  orbit  for  a  globular 
cluster,  in  order  to  arrive  at  some  typical  value  for  the  velocity  of  a  cluster 
(relative  to  the  galactic  center)  as  it  traverses  the  disk.   In  reality  the 
orbits  are  not  closed  as  we  will  see  in  Chapter  V„   However,  for  the  high  e  orbits 
with  which  we  are  concerned,  the  foregoing  simplifications  will  give  an  adequate 
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estimate.   We  Intend  to  use  this  typical  value  in  our  analysis  of  the  probable 
reaction  between  clouds  of  the  cluster  and  the  interstellar  medium  of  the  disk 
during  transit,,   We  must  recognize  that  any  such  analysis  is  by  necessity, 
order  of  magnitude  accuracy,,  We  must  realize  that  for  each  individual  transit., 
since  the  disk  clouds  rotate  with  the  circular  velocity  of  the  disk,,  the  rela- 
tive velocity  for  collision  can  vary  from  about  zero  to  approximately  600  km/sec. 
For  purposes  of  our  analysis,,  we  will  select  300  km/sec  as  the  mean  relative 

o 

velocity-   We  use  2  x  10  years  as  the  period  for  the  mean  orbit <> 
Stability  of  Clouds  and  Globules-  We  have  adopted  Roberts ' s  cloud  as  being  a 
spherical  object  of  260  solar  masses  and  2  pc  diameter-   It  is  interesting  to 
consider  the  gravitational  stability  of  such  an  object-   For  this  purpose  we 
use  Jean's  gravitational  instability  criterion  (1928),  and  find  a  rough  limit 
that  a  self  gravitating  object  must  satisfy  in  order  to  be  gravitationally 
stable  and  not  break  up  under  the  result  of  its  own  thermal  motions-   In  apply- 
ing any  such  criterion,,  we  assume  the  cloud  to  be  isolated  in  Space-   If  there 
were  a  surrounding  medium  at  a  finite  pressure,  then  the  minimum  temperature 
for  stability  would  be  larger-   The  most  convenient  form  of  this  criterion,, 
which  can  also  be  derived  from  the  virial  theorem,  is:   N>Ncr  =  2.5x10  .T/lVr 
This  signifies  that  a  body  of  mass  M  (in  solar  masses)  at  temperature  T  must 
have  a  density  of  more  than  W^cm  ^  of  hydrogen  atoms  to  be  stable-   For  the 
260  O  cloudy  we  find  the  density  of  hydrogen  atoms  to  be  about  2700  cm  .   For 

is  cloud  to  be  gravitationally  stable ,  its  temperature  would  have  to  be  below 
*+0     to  50  K-   Such  a  low  temperature  is  certainly  feasible  in  an  HI  region 
(neutral  hydrogen)  particularly  when  appreciable  obscuring  dust  is  present, 
and  there  is  a  strong  possibility  that  any  such  object  is  stable.   On  the  other 
hand,  the  0.1  O  globules  fall  short  of  the  critical  density.   However,  Idlis 
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has  shown  that  these  objects  should  have  a  lifetime  of  approximately  lCr   years 
when  considering  their  thermal  diffusion,,  so  the  globules  are  probably  long- 
lived  objects  in  their  present  state., 
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CHAPTER  III 
ENVIRONMENT  WITHIN  THE  GLOBULAR  CLUSTER 

If j  in  fact,,  we  are  convinced  that  the  clouds  and  globules  do  exist  in 
certain  globular  clusters,,  then  the  next  matter  of  interest  is  an  analysis  of 
the  environment  within  these  clusters  which  encourages  the  formation  and 
existence  of  such  objects*   Of  principle  interest  is  an  estimate  of  the  amount 

0 

radiation  from  the  stars  beyond  the  Lyman  limit  (912  A),  which  would  ionize  the 
hydrogen,,   Again ,  we  use  the  observational  data  available  for  M5« 

Star  counts  by  Sandage  for  M5  (195*1-)  have  shown  that  the  group  of  stars  brighter 
than  absolute  visual  magnitude  +3*5  (l°25  O)  is  slightly  concentrated  toward  the 
center  of  the  cluster.   It  is  in  this  group  of  more  massive  stars  above  the  main 
sequence  that  we  expect  to  find  those  actively  evolving  toward  the  white  dwarf 
stage  by  ejecting  material.  We  might 5  therefore,,  expect  to  find  this  matter  ini- 
tially with  the  same  distribution  as  these  bright  stars.   Generally^,  this  would 
mean  a  spherical  distribution^  somewhat  concentrated  toward  the  center.   More 
specifically 3   we  can  give  Roberts's  estimate  that  a  maximum  of  about  1000  solar 
masses  of  interstellar  material  is  available  between  passages  through  the  galactic 
disk.   He  makes  this  estimate  for  a  recent  epoch  by  first  finding  the  number  of 
cluster  stars  (AN)  lying  above  the  main  sequence  for  M3,  using  the  luminosity 
function  of  Sandage  (1957)°   He  also  computes  the  time  (At)  necessary  for  the  star 
to  evolve  through  this  region  of  the  color -magnitude  diagram.   Then  after  finding 
Av».j  the  average  mass  loss  suffered  by  a  star  during  evolution^  (using  1.2  0  for 
initial  mass  and  0.6  0  for  final  white  dwarf  mass)  he  estimates  the  mass  loss  of 
the  stars  of  the  cluster  for  a  half  period  of  its  galactic  orbit  (^T)  by  the 
product :    AN  .  C^-T^f)  •  ir.   Since  more  than  half  of  these  bright  stars  lie 
within  8  pc  of  the  center  of  our  M3  model,,  we  assume  as  a  first  approximation  that 
we  have  700  O  of  diffuse  material  spread  evenly  within  this  radius. 

This  approximation  allows  us  to  make  use  of  the  concept  of  "Stromgren  Spheres" 
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i  Stromgren  (1939)  to  show  the  extent  of  ionized  hydrogen  around  stars  of  types 


05  to  AO,  assuming  that  the  surrounding  material  is  distributed  uniformly  and  of 

nay  - 

-2/3 


_3 

density  1  H  atom  cm  .   These  tabulated  radii  may  then  he  modified  for  regions 


of  different  concentration  by  multiplying  by  N 

To  make  this  estimate  for  our  cluster  model,  we  use  the  computations  of 
Minkowski  and  Osterbrock  (1959)  f°r  M3-   They  found  that  about  15  hot  blue 
stars,  of  effective  spectral  type  05  and  of  considerably  fainter  absolute  magni- 
tude (M  =  +2 ..5)  than  ordinary  Population  I  stars  of  the  same  effective  temper- 
ature,  would  be  responsible  for  most  of  the  ionization.   The  uncertainty  of  the 
data  is  illustrated  by  the  fact  that  only  one  of  these  05  stars  has  been  sighted 
and  the  rest  are  inferred  from  an  incompleteness  factor  computed  by  Sandage  (195^-)° 
These  15  stars  are  distributed  within  28  parsecs  of  the  center  of  the  cluster „ 
Since  we  are  concerned  in  our  calculation  with  the  volume  within  8  pc  of  the 
center ,  we  correct  for  the  lower  number  of  stars  by  noting  that  for  the  Oort- 
van  Herk  model  (d),  the  volume  within  8  parsecs  contains  about  0,6  of  the  mass 
of  the  stars  brighter  than  M  =  +3»5  contained  within  28  pc.   So  we  correct  the 

Minkowski -Osterbrock  equation  for  the  final  total  ionized  volume  ,  V  =  ,  -^  (pc  ) 

6  ^ 

where  N  number  density,  cm"  ,  of  H  atoms.,  and  use  V  =        ( pc  ),   Stromgren 

IT 
(19^-8)  has  given  a  correction  to  his  tabulated  radii  of  ionized  spheres  for 

the  dilution  effect  for  stars  that  lie  outside  the  interstellar  medium.  Using 

this  correction,  we  find  we  can  neglect  the  contribution  of  those  blue  stars 

outside  of  8  pc  for  our  order  of  magnitude  estimate. 

Now,  using  a  total  average  mass  of  hydrogen  atoms  in  diffuse  matter  within 

8  pc  as  700  O,  we  find  the  average  concentration  of  H  atoms  is  N  -   15  cm  ,   At 

this  density,  we  find  that  the  whole  volume  within  8  pc  would  be  ionized.   This 

esult  is  for  a  uniformly  distributed  medium,  and  does  not  take  into  account  that 

any  material  would  be  concentrated  toward  the  center  as  indicated  by  the  potential 

field  in  Figure  1,   Also  there  co-old  and  probably  would  be  turbulence  and  density 
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fluctuations.  Nonetheless ,  this  particular  environment  does  not  appear  conducive 
to  the  formation  of  HI  regions,,  particularly  of  the  size  of  260  G>„   Ionized 
hydrogen  has  a  kinetic  temperature  of  approximately  10,000  K,  which  gives  a 
random  thermal  velocity  of  about  16  km/sec  that  would  cause  it  to  escape  from 
the  cluster  quickly. 

In  order  to  form  a  general  impression  of  the  effect  of  more  interstellar 
material  within  a  globular  cluster ,  we  will  hypothesize  that  3^000  O  of  mater- 
ial is  available,,  which  would  be  the  case  if  the  material  remains  undisturbed  irj 
the  cluster  for  a  time  interval "  equal-  to  three  mean  orbital,  periods*. 

About  2,000  O  of  this  would  be  within  8  pc,   We  will  again  assume  it  to  be 

1  1x10    3  3 

spread  uniformly.   Using  the  relation  V  =     a  (pc  )  with  N  now  equal  k2   cm  , 

we  find  that  less  than  one-third  of  the  volume  within  8  pc,  or  630  pc   ,  is  now 
ionized.   We  could  visualize.,  in  this  particular  case,  regions  of  ionized  mater- 
ial or  emission  nebulae  surrounding  each  hot,,  blue  star,,  separated  by  larger 
volumes  of  HI  gas.   For  the  case  between  the  two  extremes  just  discussed,,  we 
might  picture  instead  a  region  of  tenuous ,  ionized  gas  escaping  around  floating., 
smaller  volumes  of  HI, 

The  2 .,000  O  of  interstellar  hydrogen  would  be   just   ionized.,  if  distri- 
buted uniformly  within  a  volume  of  11  pc  radius.   It  is  interesting  for  this 
case  to  estimate  the  time  interval  for  ionizing  the  entire  2,000  O,  assuming 

o 

that  all  quanta  beyond  the  Lyman  limit  (912  A)  are  fully  absorbed  by  the  hydrogen 
atoms.   To  make  this  approximation,  we  must  first  calculate  the  radius  of  the 
sub-luminous  0-5  type  star,  which  is  the  principal  ionizing  agent.   We  use  a 
formula  from  Russell,  Dugan,  and  Stewart  ( 1938)  Mv  =  9^°°  _  ^  J^  r_q,08, 
om  Stromgren,  we  adopt  an  effective  temperature  for  the  05  star  of  79,000  °K. 
olving  for  R,  we  get  O.k   solar  radii. 
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The  energy  density  at  the  surface  of  a  star  of  effective  temperature.,  T, 
is  given  "by  the  Planck  law  asj 

gTThV3     l 

To  find  the  energy  density  in  the  frequency  range  "beyond  the  Lyman  limitv  ' } 
we  integrate  the  above  expression; 


u  =  8TTh 


r°  (eh^M)    V3dV 


c3 

At  these  high  frequencies ,  we  make  use  of  the  Wien  approximations 
Integrating,  we  get 


u      c, 


o  t  v.-  e- h^)  [i  «  ck)  ♦  <  (sap  *asr 

Now,  to  find  the  intensity  at  the  star's  surface,  we  use  the  relations 

LC^Kec^  =  4TTITI  L  -  1.7  -  I037  e^/iec 

This  is  the  luminosity  beyond  the  Lyman  limit  for  one  05  star  of  My  =  +2.3° 
We  can  check  these  figures  by  giving  an  approximate  B, C.  (bolometric  correction) 
for  an  Op  star.   From  Allen  (1955),  BC  =  -k2.$   +  10  log  T  +  29,000/t,  and  for 
79,000  °k^  BC  =  60 9.   Using  the  value,  K         =  k.62   for  the  sun  and 

L  :=  3„86xlO  ergs/sec,  we  find  L  ^  1.9x10   ergs/sec.   This  figure  indicates 
that  our  calculated  luminosity  1  beyond  the  Lyman  limit  is  of  the  correct  order 
Of  magnitude. 

Within  11  parsecs,  we  have  approximately  11  (05)  stars,  so  that  the  total 
energy  output  absorbed  by  2000  O  of  hydrogen  is  1.9  x  10^°  ergs/sec.   Each  atom 
of  hydrogen  absorbs  13° 5^  ev  in  order  to  ionize.   In  2000  O,  there  are  approx- 


6o 
imately  2,k   x  10   hydrogen  atoms .   This  means  we  require: 

2 A   x  10°°  x  13.5^  ev  x  1.6  x  10"12  ergs/ev  =  5»2  x  10  9  ergs 

38 
to  ionize  all  the  material.   At  the  energy  output  rate  of  1*9  x  10   ergs/sec, 

the  ionization  lifetime  is;    „  _     qq 

5.Z  «|0  _  ev-«^5 


|^<|0^e^/5ec  *3*/07^>qear    * 

Such  a  short  ionization  time  (this  is,,  of  course,  a  minimum  limit  because 
of  neglect  of  radiation  losses)  would  suggest  that  the  interstellar  material 
is  almost  immediately  ionized  as  it  appears  and  quickly  "lost  from  the  cluster,. 
Hcwever,  all  of  our-  estimates  have  been  based  on  the  use  of  Stromgren  spheres 
for  a  group  of  stars  with  the  interstellar  material  uniformly  distributed.   In 
reality,  we  would  not  expect  this  to  be  the  case.   Space  densities  derived  for 
M3  from  Sandage's  counts  for  stars  brighter  than  My =  +3°  5  show  that  these 
brighter  stars  have  a  distribution  proportional  to  R  .   King  (1962)  has 
recently  given  an  empirical  law  for  projected  densities  of  the  inner  regions  of 
globular  clusters  substantiating  this  type  of  distribution.   Assuming  that  the 
initial  space  density  of  the  material  follows  the  same  law  as  the  brighter 
stars,,  we  conclude  that  the  initial  distribution  of  matter  would  be  highly  con- 
centrated toward  the  center.   If  we  also  assume  that  the  few  stars  responsible 

-3 

for  most  of  the  ionization  follow  the  same  R  ^  law,  then  we  can  see  that  the 

central  regions  provide  the  likeliest  environment  for  breeding  HI  clouds.   The 

volume  of  ionization  varies  directly  as  the  total  luminosity  of  the  ionizing 

p 
stars  and  inversely  as  N  .   Since  the  total  luminosity  varies  directly  as  the 

number  of  ionizing  stars,  we  finally  can  write  that  the  percentage  of  ionized 

volume  varies  as  ^**/n  =  — .   In  other  words,  the  outer  regions  tend  to  ionize 

more  easily.    The  foregoing  arguments  do  not  attempt  to  provide  any  mechanism 

to  account  for  the  actual  formation  of  HI  clouds,  but  merely  offer  a  general 

analysis  as  to  the  most  favorable  regions  for  such  formation. 
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An  investigation  "by  Spitzer  (19^2)  shows  that  for  a  spherical  self- 
gravitating  system  the  interstellar  material  retained  "by  a  potential  field  ■would 
tend  to  concentrate  toward  the  center  because  of  the  field  and  the  effects 
of  inelastic  collisions.   This  effect  would  accentuate  the  initial  concentra- 
tion,.  The  fact  that  Roberts's  clouds  are  found  some  distance  from  the  center- 
seems  to  pose  an  anomaly.   This  point  will  be  elaborated  further  in  Chapter  IV. 

It  is  evident  from  the  foregoing  discussion  that  one  of  the  biggest 
detriments  to  the  formation  of  HI  regions  Is  the  presence  of  hot,,  blue  0  stars 
with  their  high  proportion  of  radiation  beyond  the  Lyman  limit.   This  would 
suggest  that  a  correlation  exists  between  the  presence  and  distribution  of 
0  stars  and  the  formation  of  HI  regions.   Any  such  correlation  would  not  be 
straightforward 9   however s   because  once  dense  HI  clouds  would  form.,  then  the 
subsequent  appearance  of  these  sub  luminous  0  stars  in  the  process  of  evolu- 
tion would  not  necessarily  ionize  dense  regions.   Nonetheless 3    it  would  seem 
logical  that  some  correlation  between  the  age  of  stage  of  evolution  of  a 
cluster  and  a  favorable  environment  for  the  production  of  HI  regions  might 
exist.   This  point  will  be  discussed  further  in  Chapter  V. 

We  next  look  briefly  at  the  kinetic  equilibrium  temperatures  we  might- 
expect  to  find  in  the  HT  clouds  of  the  cluster.   A  series  of  papers  by  Spitzer 
and  Savedoff  (19^8^  19^  1950^  1955)  have  concluded  that  the  equilibrium 

temperature  of  the  HI  clouds  in  the  galactic  plane  is  nearly  independent  of 

6  -*> 

its  density  in  the  range  from  about  1  to  10  H  atoms  cm  „   This  conclusion 

assumes  that  the  gas  has  sufficient  time  to  reach  equilibrium  with  its  asso- 
rted radiation  field. 

Before  analyzing  the  HI  regions  In  a  globular  cluster,,  we  first  might- 
discuss  some  of  the  more  important  energy  gain  and  loss  mechanisms  that  are 
involved.   These  were  discussed  by  Eddington  (1930)  and  later.,  in  more  detail 
by  Spitzer  and  Savedoff.   Throughout  the  discussion,  we  should  remember  that 
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we  are  dealing  with  a  medium  not  in  thermodynamic  equilibrium,  and  in  a  dilute 
radiation  field.   However,  the  frequent  elastic  encounters  between  the  particles 
do  tend  to  establish  an  equipartition  of  the  kinetic  energy  of  translation,  and 
the  resulting  approach  to  a  Maxwellian  velocity  distribution  can  be  defined 
by  an  effective  kinetic  temperature. 

We  can  assume  that  there  is  virtually  no  radiation  in  HI  regions  beyond 
the  Lyman  limit,,  912  A.   The  primary  heating  mechanism  is  the  superelastic 
collisions  of  electrons  with  those  atoms  having  ionization  potentials  below 
13.54  volts  and.  which  are  relatively  abundant;,  i.e.   C,  Mg,  S,  Fe,  Al. 

Heating  by  cosmic  rays,  at  least,  in  the  solar  neighborhood s    is  probably  negli- 

=3 

gible  for  HI  regions  of  1  atom  cm  '   or  denser.   The  energy  loss  mechanisms 

are  (l)  inelastic  collisions  of  electrons  with  ions  (particularly  CII  and  Sill)., 
(2)  inelastic  collisions  between  H  atoms  and  grains  (3)  inelastic  collisions 
between  electrons  and  H2  molecules,  and  H  atoms  and  H2  molecules. 

While  Spitzer  and  Savedoff  do  not  consider  time -dependent  sources.,  Kahn 
(1955)  estimates  that  collisions  between  separate  clouds  could  be  even  1 
more  effective  than  photoionization  of  the  atoms  as  a  heating  mechanism.   By 

their  analysis,  Spitzer  and  Savedoff  arrive  at  an  equilibrium  tempera  tore 

o  -3 

between  30  and  100  K  for  HI  clouds  of  density  10  cm   or  more.   At  low 

densities  between  clouds  where  radiative  cooling  processes  are  ineffective, 
temperatures  upward  of  500  °K  are  postulated.   Kahn  shows  that  cloud  temper- 
res  can  be  expected  to  drop  very  rapidly  down  to  500  °k  between  collisions, 
and.  more  slowly  below  500  K.   These  analyses  are  observationally  substantiated 
by  the  21 -cm  measurements  which  yield  a  harmonic  mean  temperature  of  about 
100  to  125  °K  for  the  HI  clouds  of  the  galactic  plane. 

In  attempting  to  apply  any  such  analysis  to  HI  clouds  in  globular  clusters, 
we  are  at  once  struck  with  the  great  uncertainties  involved.   Certainly,  we  can 
generally  assume  that  the  chemical  composition  of  the  medium  in  the  clusters  has 
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over  proportion  of  the  heavy  elements  than  does  gas  in  the  galactic  plane, 
tra  of  the  globular  clusters  are  interpreted  to  indicate  metal  to  hydrogen 
ratios  from  l/lO  to  l/lOO  that  of  the  sun,  A  lower  percentage  composition  of 
the  elements  with  low  ionization  potential  would  reduce  the  contribution  of  one 
heating  mechanism  for  HI  regions  and  perhaps  lead  to  lower  equilibrium  tempera- 
tures. 

We  can  also  expect  the  radiation  density  inside  a  globular  cluster  to  differ 
from  that  of  the  solar  neighborhood,   We  can  estimate  this  value  quite  readily,, 
We  use  the  following  notation  in  calculating  the  approximate  radiation  density 
inside  Messier  3  2 

u  =  radiation  density  (ergs/ cm )„ 

r  =  distance  from  center  of  cluster  ( cm) , 

c  =  velocity  of  light.  (3  x  10  '  cm/sec) . 

T  =  average  luminosity  (ergs/sec)  of  those  cluster  stars  brighter  than  M  =  +3°5< 

L  =  total  luminosity  of  cluster  (ergs/sec) 

*/(r)  =  density  of  stars  in  cluster  brighter  than  M  =  +3„5  ( stars/parsec  ) 
at  a  distance  r„ 
From  the  Oort-van  Herk  data  for  model  ( d)  of  M3<,  we  use  the  value  3,3  for  LM  and 
190,000  solar  mass  for  the  total  mass  (09)  of  the  cluster.   Using  the  value 
3»86xlCr   ergs/sec.,  for  the  total  luminosity  of  the  sun  and  neglecting  the  probably 
small  bolometric  correction  for  the  clustery  we  find: 

L  =  (3.86  x  1033  ergs/sec),  (190,000),  (3» 3)  =  2,i)-3xl039  ergs/sec, 

the  purpose  of  determining  the  radiation  density  of  the  cluster,  we  consider 

that  most  of  the  light  comes  from  those  stars  brighter  than  M  =  +3,5,  so  that  the 

density  distribution  of  these  stars  can  be  equated  to  the  distribution  of  light 

sources.   For  a  single  star:    u  =  --   at  any  point  a  distance  r  from  the  star, 

kitr   c 

Now,  to  account  for  the  contribution  of  radiation  of  all  the  other  stars  to 

this  single  point,,  we  integrate  over  all  the  cluster,  or  in  this  case  out  to 
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100  pc,  the  visible  limit. 

So,,  at  a  single  point  in  the  cluster,  we  have: 


\oo 


AVJi 


M-TTy^Q     V  /  /  e>  C_~ 


For  the  total  luminosity  of  the  cluster,  we  have  the  following  expression: 

(2)       L=    [J^IffrWc/r 

Dividing  (1)  by  (2),  we  get     _  _         C'^-Vclt^ 


Mttc   f  ,00-r*Vdr" 


G 


Our  final  expression  for  radiation  density  is  them 

r  \    .  ■  _  L/  f0""Vcfr 

Figures  (2)  and  (3)  from  Oort-van  Herk  depict  the  space  densities  (V)  of  the 


group  of  stars  brighter  than  M^=  +3°  5  versus  cluster  radius  (r).  These  graphs 
are  based  on  star  counts  for  Messier  3  by  Sandage  and  von  Zeipelo  In  order  to 
find  u,  we  evaluate  equation  by  carrying  out  the  integrations.   We  split  up  the 


■  erval  from  0  to  100  into  3  segments;   (l)  integrating  numerically  by  Simpson's 

r 


Rule  from  0  to  k   psc,  (2)  represent  analytically  asA/-— J !  from  k   to  7  parsecs, 

,  x    /   *» 

and  (3)  as /y  =  — ~-   from  7  to  100  parsecs. 

Upon  evaluating  the  two  integrals,  and  using  L  =  2»^3xKr^  ergs/sec,  we 
find  that  u  =  2»9  x  10 ~10  ergs/cm^. 

This  value  can  be  compared  with  the  density  of  dilute  black  body  radiation 

as  calculated  for  the  solar  neighborhood;   5° 24  x  10"  ^  ergs/cm  (Dunham,  1939) 

-13        3 

or  7« 67  x  10"   ergs/cm  (Eddington,  1930),  for  which  the  black  body  temperature 

of  a  particle  in  space  is  3«l8  °K» 


19 


k  -15     /  3   ^ 

Using  the  formula  u  =  akT  f   where  a  =  7.57x10  ^  ergs/cm  deg  and  k  is  the 

absorption  coefficient s   we  can  solve  for  the  effective  temperature  of  the  dust 

in  the  cloud  if  we  can  find  an  expression  for  ko 

We  know  from  Roberts  that  the  260  O  cloud  absorbs  approximately  5  magnitudes 

in  the  visible  light  ( Asi5xl0~^cm).   From  Wden's  displacement  law.,  A,  ->  -  =■——. 


max 


We  can  solve  for   \         for  the  effective  temperature. 

max  c- 

We  also  know  that  k  varies  as  l/\,    so  k  =  <„.  .^3,°...  ?.,,. .,  ,  . 

max 

25x"l  0  vxT  25x10    5 

We  therefore  have  k  =  ^         3    giving  us  u  =  a   >  _    T  .   Substituting 

. 2y  « 2y 

our  value  of  u  =  2.9x10"  J   ergs/cm  and  solving  for  T,  we  find  TCi  34  °K. 

The  observational  evidence  supplied  by  Roberts  (1959)  by  his  21 -cm  measure- 
ments of  WD-3  implies  an  equilibrium  temperature  of  100  °K  or  below.   Hence.,  we 
are  led  to  conclude  that  the  temperatures  of  dense  HI  clouds  in  globular  clusters 
are  at  least  as  cool  as  HI  clouds  in  the  galactic  plane. 
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CHAPTER  IV 
PASSAGE  OF  GLOBULAR  CLUSTERS  THROUGH  THE  GALACTIC  DISK 

We  now  turn  to  an  investigation  of  the  interaction  between  gas  clouds  and 
globules  in  the  cluster  with  the  medium  lying  in  the  galactic  plane ,   Previous 
investigators  have  simply  stated  that  any  gas  present  in  the  cluster  will  he 
ejected  by  each  passage  through  the  disk. 

The  galactic  disk  consists  of  a  thin  layer  of  gas  and  dust  distributed  in 
a  non-uniform  fashion.   There  are  individual  clouds  and  spiral  arms  separated 
by  regions  of  low  density..   Densities  for  each  feature  ares   typical  individual 
clouds  -  10  atoms  cm"-^  spiral  arms  -  1  atom  cm  ,  space  between  clouds  -  0,1  to 

.1  atoms  cm  '  ,  As  previously  argued.,  we  assume  most  of  the  hydrogen  to  be 
neutral  for  purposes  of  the  collision  process.   The  ionized  gas  in  the  disk  is 
concentrated  around  the  individual  0  and  B  stars. 

Since  we  are  interested  in  the  processes  that  tend  to  act  on  the  neutral 
hydrogen  clouds  and  globules  in  the  clusters ,  we  can  see  that  the  collision  to 
be  studied  is  the  HI -HI  impact. 

As  long  as  the  mean  free  path  for  particle  collision  is  small  compared  to 
the  dimensions  of  the  cloud  or  globule,  we  can  treat  the  gross  dynamics  of  the 
encounter  by  evaluating  the  momentum  transfer  to  determine  if  the  clouds  and 
globules  of  the  cluster  are  ejected  from  inertial  force  considerations.   For 

o 

collisions  of  neutral  H  atoms,  we  can  use  the  geometrical  radius  (~10  cm)  to 
compute  the  collision  cross-section,  giving  (J=  3x10"  cm  (Parker,  1958),   For 


3-3  3-3 

-Jrcm  ^1    for  the  clouds,  W^1^™  ' 
H  H 

=6 


the  globules,  N — 109cm  |  for  the  clouds,  NorlCTcm  .   Using  the  relation, 


Mmean  free  path)-  i/no%  we  find  X   glob=  10"  pc  and  X   cloud— 10 "  pcj  with  both 
values  considerably  smaller  than  the  dimensions  of  the  object  involved, 

Kahn  (1955)  Has  considered  the  statistics  of  size  and  distribution  of  clouds 
in  the  galactic  plane  by  using  the  observational  evidence  that;   (l)  a  line  of 
sight  of  1  kiloparsec  intersects  an  average  of  10  clouds,  and  (2)  clouds  occupy 
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about  5$  of  the  volume  of  the  galactic  disk.,   He  concludes  that  the  average 

-3 
cloud  is  3-75  parsecs  radius^  with  a  density  of  20  H  atoms  cm   and  that  the 

number  density  of  clouds  in  space  is  2.%>xl0"  pc"  .   This  result;,  of  course }    is 

characteristic  of  the  solar  neighborhood,,   We  readily  see  that  the  mean  free 

path  for  collision  is  about  100  parsecs . 

From  21  cm.  observations.,  Rougoor  and  Oort  (i960)  have  described  the  general 
galactic  disk  between  5  and  10  kpc  from  the  galactic  center  as  being  220  pc  thi  .-K 
with  a  mean  density  of  1  H  atom  ci"-5,   The  disk  inside  3  kpc  seems  to  be  much 
thinner s   about  120  pe^,  and  the  mean  density  is  giver,  as  0,!j  H  atoms  cm  .   They 
comment  that  the  region  between  3  kpP  and  the  central  disk  at  600  pc  appears  to 
be  virtually  empty. 

From  our  discussion  in  Chapters  II  and  III.,  we  note  the  following  informa- 
tion concerning  the  clouds  and  globules ;   (l)  In  general 5  we  would  expect  to 
find  HI  regions  toward  the  center  of  the  cluster  in  an  isolated  cluster. 
(2)  From,  the  incomplete  information  now   available 3    clouds  seem  to  appear  at  the 
edge  of  the  clusters  (minimum  of  8  to  10  pc)  and  globules  toward  the  center  of 
a  cluster  (in  the  one  example 5  M4-). 

Let  us  assume  that  a  globular  cluster  contains  a  260  O  cloud  of  Roberts's 
description  at  the  center  of  the  cluster.   Let  us  then  analyze  the  possible 
effects  of  a  collision  with  the  galactic  disk.   In  this  analysis  we  will  use 
two  different  treatments  for  the  two  different  types  of  impact  involved.   If  we 
assume  that  the  cluster  traverses  a  mean  distance  through  the  disk  of  300  parsecs., 
then  at  300  km/sec  (where  1  km/sec  Cil  parsec/lO  years)  we  can  conclude  the 
time  of  transit  is  about  10  years.   However.,  a  cluster  cloud  would  transit  one 
of  Kahn's  typical.,  dense  clouds  in  the  galactic  plane  (diameters 7. 5  pc)  in 
about  10  years.   If  we  compare  these  two  characteristic  times  with  the  approxi- 
mate orbital  period  of  an  object  in  the  potential  field  of  the  cluster.,  10'  years,, 
we  see  the  need  for  different  treatments. 
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gientum  Transfer 

For  the  case  of  impact  with  a  discrete  cloud  in  the  disk,  we  consider  an 
impulse  to  he  acting  on  the  cluster  cloud.   We  will  also  assume  that  the  colli- 
sion is  completely  inelastic,  the  colliding  portion  of  the  disk  cloud  remaining 
with  the  cluster  cloud. 

For  the  case  of  impact  with  the  intercloud  medium  (0.1  to  0.01  H  atoms  cm  } 
of  the  disk,  we  compute  the  momentum  transfer  per  unit  time  of  the  colliding 
medium,  which  is  the  force  acting  on  the  cloud  or  glohule  of  the  cluster. 

We  will  see  that  one  effect  of  the  collision  with  a  disk  cloud  can  be  to 
knack  a  cluster  cloud  out  of  the  cluster  center  and  into  a  rectilinear  orbit. 
On  the  other  hand,  the  force  exerted  by  the  medium  on  a  cluster  cloud  effectively 
shifts  the  whole  orbit  in  the  direction  away  from  the  direction  of  motion. 

The  simplest  manner  in  which  to  visualize  this  problem  is  to  analyze  it 
with  respect  to  the  potential  well  of  the  cluster  and  to  consider  the  motion 
as  a  one  dimensional  oscillator.   We  use  the  following  notation: 

V  =  velocity  of  globular  cluster  relative  to  the  Galaxy  (300  km/sec )„ 

R  =  distance  of  cloud  or  globule  from  center  of  cluster. 

r  =  radius  of  cloud  or  globule  of  cluster. 

L  -  diameter  of  disk  cloud  ('7-5  Pc) 

n  =  number  density  of  H  atoms  in  galactic  disk  (cm  ). 

-2k 
m =  mass  of  hydrogen  atom  (1.67x10   gm). 
n 

(T/|  =  mass  of  cloud  in  cluster  (260  O). 

$(R)  =  potential  as  a  function  of  R  in  the  cluster  (km/sec)  . 
A$  =  potential  gradient. 

N  =  number  density  of  H  atoms  in  globular  cluster  cloud. 
F  -   Momentum  per  unit  time  or  force  applied  to  cloud  in  cluster  resulting 

from  transit  of  intercloud  medium  of  disk. 
E  =  Total  energy  per  unit  mass  of  cloud  in  globular  cluster. (km/sec)  „  (with 

respect  to  cluster  frame  of  reference). 
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Figure  ( l)  depicts  our  detailed  example  of  the  globular  cluster  -with  a  260  O 
cloud  at  the  center  of  the  cluster  on  its  first  impact  with  the  disk. 

First,  we  derive  a  general  formula  for  computing  the  force  per  unit  mass 

resulting  from  momentum  transfer  per  unit  time  from  the  intercloud  medium,  to 

P  P 
the  cluster  cloud  on  transit  of  the  disk.  We  have  that;   F  =  itr  .V  .h.il,(  dynes) 

3 

nn  kitr 


N.  il(  grams ) 

L   (  ^m/RPf'  ^ 

4  r.N 


(1)  F/f^JI-^  (em/sec2) 


2  ""   3 


5  =3  / 

For  our-  particular  case,  where  r  =  1  pc  and  N  =  2.7-lO^cm  ,  V  =  300  km/sec  and 

n   =   0.1  crrT5_,  equation  (l)  gives  us;   F./p?  =  8. 2x10 °9 cm/sec2  or  2.5[— ■)  pc"1, 
We  can  also  give  for  our  example,  the  expression  for  the  force  as; 

(2)  F(  dynes)  =  lj-.25xl028.n(om~5) 

or       F  =  4.25xl02^  dynes  (for  n  =  O.l). 

We  have  the  following  relationships  between  total  energy  per  unit  mass  (E),  work 

F.CRp-Rj 

done  on  the  cluster  cloud  per  unit  mass  potential  energy  per  unit. 

mass  ($)  and  kinetic  energy  per  unit  mass  {2.^) ,  where  t\f  =   velocity  of  cloud  in 
its  orbit  with  respect  to  the  globular  cluster; 

(3)  E  =  E_  +x?/R^-R1)  where  E  is  the  total  energy  in  an  orbit  at  time  ts 

2     1  \  *7  2   1         1 

and  E  is  the  total  energy  at  some  time  t  later.   As  shown  in  figures  (l)  and  (2) 
R  values  in  the  direction  away  from  the  direction  of  impact  are  positive.   We 
also  have  the  following  identity  relating  our-  energies; 

E  =  O+^T  (where  Vis  the  velocity  of  the  cloud  in  its  orbit  relative  to 

the  frame  of  reference  of  the  cluster) 

2 
We  can  see  that  the  value  of  YJ(f\ «  — ~     determines  to  what  extent  an  orbit  is 

displaced  to  the  rear  because  of  the  steady  force  exerted  by  the  intercloud 

6 
medium  during  about  10  years.   In  figure  (2)  we  have  depicted  the  effect  of 

force  per  unit  mass  exerted  by  both  a  0,1  cm"^  medium  (F/(f)  =  2. 5( km/sec)  pc"  ) 

and  a  0,01  cm"^  medium  (f/(T)=  0.25).   We  depict  the  extreme  case  for  both  these 
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media  In  order  for  a  cluster  cloud  (of  Roberts's  description)  to  be  retained 

-3 

in  the  cluster  for  an  indefinite  period  of  time,,   For  the  0.1  cm  medium.,  the 

limiting  orbit  is  from  -k   pc  to  +1^  pcj  for  n  -   0-01  cm~-^  it  is  from  -20  pc 
to  +50  pc. 

We  should  note.,  however .,  that  for  our  particular  case  where  the  time  of 
transit  of  the  disk  is  10  years  and  the  orbital  period  in  the  cluster  is  ~10' 
years,  that  the  chance  of  a  cloud  escaping  due  to  the  reaction  of  the  medium  alone 
is  extremely  remote.   Once  the  cluster  emerges  from  its  transit  of  the  disk., 
then  the  F/jf]  slope  is  practically  zero  and^  as  long  as  the  cloud  has  a  total 
energy  less  than  138  (km/sec)  9    it  will  not  escape. 

However ^  the  reaction  caused  by  collisions  with  discrete  clouds  in  the 
galactic  disk  of  the  order  of  n  =  10  cm  ^  is  extremely  significant  for  our 
discussion.   To  calculate  the  energy  change  per  unit  mass  (AE)5  we  use  the 
principle  of  conservation  of  momentum.   For  the  general  case  we  have  for  the 
s  of  material  intercepted  by  the  cluster  cloud; 


(Y|.=  L.«.^ 


rui^ 


lffiH 


The  mass  (00)  of  the  cluster  cloud  is; 

V  Is  the  relative  velocity  of  impact  (300  km/sec). 

A^is  the  velocity  change  imparted  to  the  cluster  cloud  relative  to  the  cluster 

frame  of  reference. 

We  then  have;(V}.av  =  (09  +P9j) .  A  W 

\/'OQ.  ' 

This  gives  the  expression  A'U  ^  .  . v  -J ,  c 

(4)   Aor= 


Now  we  derive  a  general  expression  for  the  energy  transfer  per  unit  mass  per 
collision.   We  haves 
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E.  =  $+2 Oft'     (where  *V.  is  the  orbital  velocity  of  the  cloud  at  time  t.) 
i       v   V       1  J  i' 

E  =  $+■§(  AT. +A AT)    (E  =  energy  after  impulse) 

This  gives  us  the  general  expression  for  energy  transfer; 

o 

(5)   AE  =  E  -E.  =  1TAU+  LA^L 
f  l    l       2 

Substituting  the  appropriate  values  into  equation  (k)   for  the  260  O  cluster 
cloud  colliding  with  the  n  =  10  cm   disk  cloud, we  find  Afl/s:6  km/sec.   From 
equation  (5)  then.,  we  have  AE  =  6  (1T.+3). 

NOW;,  with  the  use  of  figure  (l)j  let  us  follow  through  in  some  detail  the 
reactions  of  a  260  O  cloudy  initially  at  rest  in  the  center  of  a  globular  cluster 5 

as  it  passes  through  a  typical  area  of  the  galactic  disk.   We  adopt  0. 1  H 

-3 
atoms  cm  '  as  the  concentration  for  the  intercloud  medium.   We  assume  that  the 

average  cluster  passes  through  $00  pc  of  the  disk.   The  mean  free  path  for 

-"5 

collision  with  one  of  the  disk  clouds  (L  =  7.5  pc5  n  =  10  cm  )  is  100  pc. 

Upon  gradual  entry  of  the  cluster  into  the  0.1  cm  v  intercloud  medium,,  we 
expect  the  cloud  to  "be  shifted  about  0. 2  pc  to  the  rear  of  the  cluster  center^ 
or  to  R  =  +0.2  pc.   The  cluster  cloud  is  likely  to  be  near  this  point  on  the 
diagram  upon  its  first  collision  with  a  disk  cloud.   Using  our  previously  calcu- 
lated value  of  A^  6  km/sec^  we  find  from  equation  (5)^  where  1Tj,=  0,,  AE  =  lfir—  J 
We  indicate  this  change  in  energy  on  figure  (l)  by  a  dotted  vertical  line  to 
point  (R   e  )  or  (+0.2;,  19 )•   At  this  point  we  assume  that  about  one  third  of 
the  transit  period  (^xlO-^years)  has  elapsed.   Through  this  point  we  draw  a  line 
with  the  proper  slope  to  represent  F/JqC^2.5  (—  )  pc"  ,   as  previously  calculated 
for  the  force  per  unit  mass  exerted  by  the  0.1  cm~^  medium  on  the  260  O  cloud. 
This  line  represents  the  locus  of  (R,E)  points  for  the  rectilinear  orbit  which 
the  cluster  cloud  now  follows. 
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Figure  1  -  POTENTIAL  WELL  FOR  M5  SHOWING  TYPICAL  PASSAGE  THROUGH  DISK 
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Figure  2    -  POTENTIAL  WELL  FOR  M5   to  100  pc 
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The  increment  of  energy s   AE  =  18  (-!——)  3   received  by  the  cluster  in  the 

S6C 

first  collision  has  given  the  cloud  an  impetus  toward  positive  R  values.  We 

now  estimate  the  distance  traversed  along  this  orbit  in  a  time  interval,, 

5 
3x10  years 5  until  the  next  collision.  We  integrate  numerically  along  the 

path  from  R  -  0.2  pc  to  an  R  value  that  corresponds  to  a  time  interval  of 


xlOv  years.  We  use  distance  steps  of  0.2  pc; 
R1+n(0.2)   ^ 


a 


j- 

We  then  find  that  the  cloud  will  be  at  (+1.6,  22)  at  the  time  of  second  collision. 

As  depicted  in  figures  (l)  arid  (2),,  we  follow  the  cloud  through  two  more 
collisions.   After  the  second  collision  and  a  time  interval  of  3x10^  years^ 
we   find  that  a  AE  =  58  (km/ see)   for  the  third  collision  raises  the  cloud  to 

very  high  energy  total  (E  —  125  (km/see)^  which  means  the  cloud  probably 
escapes. 

We  note  from  figure  (l)  that  if  the  cloud  does  not  undergo  the  third 
collision.,  it  will  find  itself  in  an  orbit 3  with  the  distance  limited  to  about 
7  parsecs.   From  figure  (2) 9   we  see  that  a  cloud  after  the  third  collision  is 
supposedly  left  in  an  orbit  extending-  out  to  about  60  to  JG   parsecs.   However ^ 

it  is  easy  to  see  that;,  even  for  a  very  small  fractional  force  which  could 

-2      -K        -3 
reasonably  be  expected  from  the  hale  medium  (10~  to  10  ~  cm  ).,  the  cloud 

would  be  lost  from  such  a  high  energy  orbit. 

While  we  have  worked  out  only  this  one  example  in  any  detail,  it  is  easy 
to  judge  the  possible  reactions  in  other  situations. 

For  instance.?  a  subsequent  pass  by  our  cluster  in  the  example  above -which 

o 

has  had  two  collisions  on  the  previous  pass  through  the  disk  10  years  before.- 
might  find  its  cloud  in  an  orbit  at  an  energy  level  of  about  E  =  70°   The 
initial  velocity 5  for  this  ease  might  vary  anywhere  from  about  t\f„    =  +12  km/see 
through  0  and  back  to  -12.   Minus  values  for  velocity  indicate  a  direction  for 
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.cloud  in  the  same  direction  as  the  motion  of  the  cluster  orbit.   Using 

xrmulae  (k)   and  (5)*  we  find  that  AE  can  have  the  following  values  for  each 

value  of  ^T„«   for  +12  km/sec,  AE  =  +90^  0  km/sec.,  AE  =  +l8j  -12  km/sec^  AE  =  -^h, 

Hence^  we  see  that  in  these  3  cases;   the  cloud  escapes^  has  its  energy  raised 

about  18  (km/sec)  }   and  has  its  energy  lowered  ^>k   (km/sec)   respect ively, 

1  .vcistically,,  we  would  expect  the  cloud  to  he  near  an  extreme  point  of  its 

orbit  and  thus  find  that  the  change  in  energy  for  an  impulse  is  small. 

Two  other  factors  are  worth  "brief  mention;   (l)  we  ^could  expect  some 

elliptic  orbits  to  ensue.for  the  reason  that  a  cluster  on  a  subsequent  pass 

could  collide  with  the  disk  at  an  angle  to  the  direction  of  motion  of  its  rectilinear 

orbits  and  (2)  we  might  expect  some  damping  of  the  orbit  during  the  period  the 

cluster  is  away  from  the  disk.,  because  of  accumulation  of  additional  medium 

from  evolution  of  stars  and  the  effects  of  inelastic  collisions. 

To  estimate  the  chance  of  a  particular  cluster  cloud  getting  through  the 

ik  Without  escaping  from  the  cluster s   we  use  the  results  from  the  previous 

ocussion  and  apply  the  Poisson  frequency  function; 

k  -a 

where  a  is  the  expectation  and  f(k)  is  the  probability  that  x  -   k.   For  our 

problem.,  we  take  a  equal  to  3  (we  can  expect  about  3  collisions  per  pass). 

We  then  evaluate  the  individual  probabilities  for  k  =  0^,  1.,  2  and  3  collisions 

transit.   We  sum  up  these  individual  probabilities 3   using  \„  f(k)  for  k  =  3.? 

find  that  the  probability  of  a  cluster  cloud  surviving  one  disk  passage  is 
aboui     Of  course  this  figure  is  for  the  average  assumed  conditions  of  the 

:  and  would  vary  for  a  different  structure  in  any  region. 
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We  should  note  that  the  foregoing  development  for  momentum  transfer 
applies  equally  well  for  the  Idlis  cloud.   That  this  is  true  is  readily 
understood  from  equations (l)  and  (k)   which  show  that  the  forces  depend 
on  the  product  rN  which  are  approximately  equal  for  the  two  cases. 

It  is  clear  that  any  original  gases  in  the  globular  cluster  are 

10 
swept  out  at  some  time  during  passages  through  the  disk  during  10   years. 

We  previously  mentioned  one  other  possible  source  of  material  for  the  cluster s 

that  of  sweeping  up  the  gases  from  the  halo.   Spitzer  gives  a  value  of 

5x10 =  H  atoms  cm""  for  the  halo.   Shklovsky  and  Pikelner  (1958)  give 

J-^1  cm  .   If ^  for  example.,,  a  Roberts  cloud  swept  a  path  during  its  orbit 

of  approximately  30  kpc  lengthy  and  retained  all  swept  matter^  the  total 

mass  accumulated  would  be  only  about  20  GX,  even  by  using  W  =  10"  „   We 

can  conclude  that  a  "snow  plow"  accumulation  is  a  relatively  small  source 

of  material. 
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; gy  Transfer 

In  addition  to  the  effects  of  moment-am  transfer,  we  must  also  consider  the 
beating  effects  of  a  collision  between  clouds  at  these  high  energies .   From  our 
discussion  in  Chapter  III  of  heating  and  cooling  processes  in  cosmic  clouds,  we 
learned  that  if  an  HI  region  becomes  ionized,  it  will  have  an  equilibrium  temper = 
ature  of  about  10,000  k  and  the  root  mean  square  velocity  of  the  protons  will 
be  about  16  km/sec.   The  Oort-van  Herk  dynamical  model  for  MJ  gives  a  velocity 
of  escape  of  16.6  km/sec  at  the  center  ranging  down  to  11.4  km/sec  at  8  pc  and 
so  on,  as  can  be  surmised  from  figure  2.   We  conclude  that  when  hydrogen  is 
raised  to  a  kinetic  temperature  of  10,000  K,  it  will  quickly  expand  and  escape 
from  the  cluster. 

Bates  and  Griff ing  (1953)  have  calculated  the  ionization  cross  section  for 
the  reaction  H(lS)+H(lS)  =  H( IS ) + l++e .   Extrapolating  their  data,  we  estimate 
this  cross  section  as  10"'"  cm  for  a  velocity  of  about  300  km/sec.   Comparing 

the  ionization  cross  section  (~10~'"  cm  )  with  the  collision  cross  section 

/    -16  2 \ 

V~10   cm  ),  we  see  that  an  individual,  hydrogen  atom  has  a  very  low  probability 

(~L0   )  of  ionization  per  collision.   The  incoming  energetic  hydrogen  atoms, 

relative  to  the  frame  of  reference  of  the  dense  (~Kr  to  10  cm  )  cloud  or 

globule  of  the  cluster,  will  be  quickly  damped  to  energies  below  13° 6  ev  as 

a  result  of  both  momentum  transfer  and  radiation  losses  by  excited  atoms.   As 

we  will  see  later,  the  ionization  will  be  mostly  confined  to  a  narrow  interface 

between  two  colliding  clouds. 

We  should  keep  in  mind  that  collisional  ionization  is  not  a  very  efficient 

process.   Alfven  (i960)  has  proposed  that  collisional  ionization  between  a 

plasma  and  a  neutral  gas  is  caused  by  electrons  being  accelerated  to  sufficient 

ergies  by  irregularities  in  the  charge  distribution  created  in  the  plasma. 

In  the  discussion  following  this  paper,  Spitzer  remarks  that  when  ionization  is 

produced  by  energetic  electrons,  the  cross  section  for  excitation  is  much  greater 

than  that  for  ionization  and  k   to  5  times  as  much  energy  goes  into  light  than 
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into  ionization.   Naturally.,  for  any  ionization  mechanism,  a  certain  proportion 
of  the  input  energy  will  be  dissipated  as  light  and  infra-red  radiation  from 
excited  atoms  and  from  grains* 

We  will  analyze  in  some  detail  the  case  of  collision  between  a  Roberts 
cloud  (N~10  cm  )  and  a  typical  cloud  in  the  disk  (n~10  cm  ,  diameter  [l]~7.5P ■")« 

First.,  we  should  note  that  at  relative  velocities  of  about  300  km/see,  the 
interacting  particles  exceed  the  velocity  of  sound  and  we  can  expect  a  compressed 
shock  front  between  two  colliding  clouds.   The  velocity  of  sound  is  given  by 
a  =\jtRT*  where  ^  5/3,  R  =  3. 3x10  \   For  T  =  100  K,  a  -  1.2  km/secj  for 
T  =  10,000  °YL3   a  =  17  km/see. 

We  now  calculate  the  total  energy  per  second  transferred  to  the  more  dense 
Roberts  cloud  as  it  passes  through  the  larger  cloud  of  the  disk.  We  have  for 
energy/unit  times 

E/unit  time  =  J( jtr  xV)x( m .  rm)xV  =  6.4x10  ergs/sec 
(where  m  =  10  cm"''  for  disk  cloud,  r  =  1  pc  for  Roberts  cloud,  V  =  300  km/sec, 

c  2  f 

-24    \  6. 4x1'.')  -2    -1 

rtn    =-.   1. 67x10   '  gm).   This  means  that  ° """ %    or  0.23  ergs  cm  sec  ""  is  the 

jtr 
amount  of  energy  input  per  unit  time  per  unit  area  which  must  be  dissipated  at 

the  interface  of  the  Roberts  cloud.   We  note  that  this  value  is  comparable  to 

the  average  stellar  radiation  per  unit  area  in  a  cluster  \ 

uc   2.9xl0~10x3xl010   _  on        -2   -1 

t—  =  — - —  1      — —  =  0.  22  ergs  cm  sec 

We  should  investigate  the  possible  cooling  mechanisms  which  are  able  to 
dissipate  this  energy  and  maintain  the  layer  at  a  reasonable  equilibrium  temper- 
ature , 

Spitzer  (1949)  has  given  the  following  formula  for  computing  the  rate  of 
energy  loss  caused  by  the  inelastic  collisions  between  hydrogen  atoms  and  dust 
grains : 

L   (  ergs/cnrsec)  =  9.4x10    n^g)N(H)T2( T-Tg) 

llg, 
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We  take  the  foil  awing  values  for: 

n(g)  -  concentration  of  dust  particles  =  10"  cm 

(where  radius  of  particles  C/LO  cm  and  density  crl  gm/cm  ) 

N(H)  -  value  for  Roberts  cloud/vlCTcm"^ 

(the  above  two  values  give  a  mass  ratio  of  dust  to  gas  of  approximately  ifo) 

o 

T  -  gas  temperature  assumed  to  be  10^000  Ko 

T   -  temperature  of  the  grains  is  assumed  very  small  (10-50  K). 

-?1        -■'5    —3 
Inserting  these  values  in  Spitzer's  formula  j,  we  find  L   =  9°^xl0~  ergs  cm  £ 

--& 

We  have  an  energy  input  by  collision  of  0,2/5  ergs  cm~2sec~  s   which  even  if  dis- 
tributed throughout  a  layer  of  several  mean  free  paths  for  collision 
(A  =  1/uW  =  1/10  „  3x10"   =  3x10  cm)  would,  be  many  orders  of  magnitude  greater 
than  the  cooling  rate. 

We  assume  that  the  impact  layer  between  the  two  clouds  is  completely  ionized 
and  at  a  temperature  of  105000  K.   Accordingly 3   we  will  have  a  recombination 
of  electrons  and  protons  with  the  accompanying  radiation  in  the  visible  and. 
infra-red.  as  the  electrons  cascade  down  to  the  ground  state.   We  use  the  follow- 
ing formula  to  obtain  an  estimate  of  the  cooling  rate  provided  by  this  mechanisms 

(&)h         (ergs  cnf-sec"1)  =  N  N  OT  ^ ( h V  )  .D 
rec.  p  e  Cn"  n 

where  N  =  proton  concentration  in  the  impact  layer  (cm  ). 
P 

N   =  electron  concentration  (cm  ) „ 
e 

t\T   =  velocity  of  electron  at  10^000  °K„ 

th 

Qm  =   capture  cross  section  to  n   state  for  hydrogen, 

(hv)   =  average  energy  radiated  per  capture  to  state  ns   and  in  cascade 

to  ground  state. 
D   =  thickness  of  ionized  impact  layer  (cm). 

From  page  91  of  .Allen  (1955).?  "we  take  the  values  given  for  the  recombina- 

th 

tion  cross  section  to  the  n  "  state  for  hydrogen  and  multiply  each  value  by 

(hV)   and  then  sum  all  products  to  gives  £<T('hV)  =  5 . 7x10 "^  eras  cm~d 

n  „22  2     n    n  D 

(where  eachfl*   is  of  the  order  of  10  cm  ). 

v  n  ' 
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For  the  velocity  of  the  electron  at  10,, 000  K.,  we  use  (IT=  6.  taclO  em/sec.   We 

consider  that  N  Osijj  and  set  the  equation  for  L    equal  to  the  energy  input „ 
p   e  re  c 

We  then  solve  for  HTDi 

KFD(V)(^  (T  (hv  )   =  energy  input  from  collision 

A  (6o'4xl07)(5„7xl0~52)  --,  0.23  ergs  em^sec"1 
fy)or  ITD  =  6o3xl022cm"5 

Having  obtained  an  expression  for  N  D^  we  now  find  N  from  considerations 
of  pressure  equilibrium  in  the  collision  process.   We  note  that  the  collision 
picture  can  be  described  in  the  frame  of  reference  in  -which  the  ionized  shock 
layer  between  the  two  clouds  is  stationary  so  that  a  disk  cloud  (n  =  10  cm™-"1) 
moves  with  a  velocity  of  about  v  =  300  km/sec.  We  can  readily  compute  the  rate 
of  transfer  of  momentum  per  unit  area  at  the  front  by  the  expression; 

n.m  ,  v2-  10  (  cm"3 ) .  1, 67. 1024(  ^n.) .  ( 3xl07)2cm2see~2 


=  1.5.1' 


~8      2   -2 
gm  cm  see- 


Since  this  value  is  the  force  per  unit  area  or  pressure  exerted  on  the 

layer  from  the  direction  opposed  to  the  direction  of  motion.,  we  can  equate  it 

to  our  expression  for  the  pressure  inside  the  layer ? 

P  =  N  kT  (where  N  =  number  of  protons  per  cm  s   k  =  Boltzmann' s  constant 
£  P 

( 1.38x10"  ergs/degree)^  and  T— 10^000  K).  We  use  N  as  the  concentration  in 

this  expression  since  the  electrons  contribute  only  about  1/V3  of  the  effective 

h     -*> 

pressure.  We  find  N  ciio  cm  . 

P 
Using  our  result  from  equation  (7)   that  N  D  =  6.3x10  cm  3   we  then  solve 

for  the  thickness  of  the  shock  layer^  finding  D  =  6.3x10  cm.   We  note  that  D 

is  larger  than  the  mean  free  path  for  collision  of  H  atoms  (\   =  3x10  cm)  in 

1    dense  cloud  and  comparable  to  the  mean  free  path  (A.  =3x10  cm)  in  the  disk 

cloud.   We  also  note  that  it  is  extremely  thin  compared  to  the  diameter  of  the 

cluster  cloud  (2  pe). 
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It  follows  that  we  have  a  situation  in  which  the  two  clouds  and  their 
interface  are  in  pressure  equilibrium,,   Using  the  expression  for  momentum 

O  O  T  T  (   2  f\  C^\   1    n  *"\ 

transfer  per  unit  area,  n.m  J",  we  see  thatsOr^      =  ~—  =  ^  ";  ° 

h  w  cluster    N      n  _3 

10 

or  0/",   .   —30  km/sec  (relative  to  frame  of  reference  where  shock  layer 
is  stationary). 

We  then  see  that  the  distance  of  travel  by  the  shock  front  relative  to 

the  cluster  cloud  during  the  time  of  collision  is  equal  to  _2_  x  L  (where 

300 

1  =  7"  5  pc)»   This  distance  of  travel  by  the  shock  front  (■'■-Co  75  Pc)  is  less 
than  the  diameter  of  the  cluster  cloud  (2„Q  pc).  as  we  would  expect. 

The  effect  of  this  impact  layer  is  to  convert  high  energy  particles  into 
radiation  in  the  visible  and  infra-red  region  at  the  spectrum.   This.,  in  turn, 
serves  to  increase  the  radiation  density  at  these  wavelengths  inside  the  cloud 
of  the  cluster.   In  order  to  determine  if  the  cloud  can  handle  this  increase  in 
radiation  density  at  these  wavelengths ,  we  analyze  the  cooling  capacity  of  just 
one  mechanism  -  that  of  the  dust  grains.   If  we  should  find  that  an  increase  in 
radiation  density  would  require  such  a  high  temperature  for  the  dust  particles 
that  they  could  not  remain  in  the  solid  state s   then  our  cloud  could  possibly 
be  raised  to  a  high  equilibrium  temperature.   We  pur-sue  this  analysis  in  a  way 
which  seeks  merely  to  set  a  limit  as  to  the  amount  of  cooling  the  dust  can  handle 
with  respect  to  the  increased  energy  input. 

From.  Roberts's  data  we  have  that  the  cloud  absorbs  about  5  magnitudes  in 

the  visible  region  of  the  spectrum  ( A^xlO^cm).   If  we  assume  that  the  dust 

a  temperature  of  about  50  °K  (at  50  °K,  the  cloud  will  be  self -gravitating), 

n.  from  Wien's  laws  A.    (cm)  =   ,   " 

max.^  "    T(oK) 

a*    (at  50  OK')^5.8xl0"5cm 
max  x        ' 

Since  the  absorption  law  for  small  (~10~5cm)  dust  particles  shows  that  absorp- 
tion varies  inversely  as  wavelength.,  we  conclude  that  effective  absorption 

-  36  - 


o                      i^xl-0       \ 
at  50  "K  Is  approximately  5  magi  ■<< *■"■  ■  j«SSLO.O^  magnitude.   This  means  that 

the  cloud  has  an  absorption  coefficient 5  (k)^5$.>  at  the  relevant  wave- 
length. \ 

max„ 

Therefore  _j  we  apply  the  Stefan -Bolt zmann  radiation  law  (with  the  additional 
factor  k)  to  our  cloudy  and  thereby  determine  how  much  energy  per  second  the 
dust  in  the  cloud  can  radiate  at  50  °Ks 

L( ergs/sec)  =  ktcc   <TT  »k 
where  r  =  1  pc^  <F  =  5. 7x10 "* erg  cm"  deg"  sec"  , 

T  =  50  °K  and  k  =  0.05. 
Substituting  these  value  S;,  we  find,  that  L»  2xl0^ergs/see  (coding  rate  that 
can  he  handled  by  a  Robert  :-  .1  >ud  with  dust,  particles  at  a  temperature  of   50  °K). 
Now  if  we  compare  this  with  the  maximum  possible  energy  input  from  the  collision 
process^  6.4x10'   ergs/sec5  we  conclude  that,  the  :leud  can  absorb  this  addi- 
tional energy  input  without  vaporizing  the  du  it. 

By  comparing  with  the  above^the  energy  input  expected  in  the  collision  of 
luster  cloud  with  the  intereloud  medium,  of  0-1  zm"J   concentration^  we  con- 
elude  that  the  Roberts  cloud  will  not  be  heated  appreciably  except  for  a  thin 
layer  at  the  collision  juncture . 

We  also  find  that  the  Idlis  cloud  can  withstand  the  thermal  effects  of 
cloud  impact  in  the  galactic  disk. 

By  a  similar  analysis  as  used  for  the  260  O  cloud.,  we  find  that  the 

k     -3  Ik 

proton  concentration  is  also  10  cm   in  the  Idlis  cloud.,  making  D^oxlO  cm. 

Using  the  expression  for  transfer  of  momentum  per  unit  area^  we  find  that  the 

ionized  layer  travels  at  0*3  km/sec  through  the  globule  (t\f  ^  i^WJ  °1U  ^ 

globule     5 

Since  this  is  below  the  speed  of  sound  (a  =  1,2  km/sec  for  the  material) ? 
do  not  have  a  shock  phenomenon  in  this  case. 
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CHAPTER  V 
DISCUSSION  OF  RESULTS  MD  CONCLUSIONS 

From  the  calculations  in  Chapter  IVS   we  can  conclude  that  a  globular  cluster 
will  not  necessarily  lose  its  interstellar  material  in  passing  through  the  galaet  i 
disk. 

We  have  seen  that.,  for  a  cloud  of  the  proportions  described  by  Roberts  and 
also  for  the  smaller 5  denser  globule  of  Idlis^  momentum  transfer  will  be  sufficient 
in  about  half  the  cases  to  eject  it  from  the  cluster 0  We  also  learned  that  a 
shock  front  forms  during  collision  with  clouds  of  the  galactic  disk  and.  that  the 
recombination  of  hydrogen^ and  radiation  in  the  infra=red  by  dust  grains  is  a  suf- 

Lent  cooling  process  to  dissipate  the  energy  of  collision  which 5  therefore,, 
will  not  cause  heating  and  expansion  of  the  cloudo 

From  considerations  discussed  in  Chapter  III,  we  concluded  that  the  center 
of  the  cluster  was  the  most  satisfactory  environment  for  HI  regions » 

Accordingly,,  it  would  seem  plausible  that  diffuse  material  accumulates  in 
.  center  of  the  cluster  and  is  excited  into  orbits  at  various  distances  on 
■sages  through  the  galactic  disk,,  We  would  not  expect  to  see  those  clouds  in 
the  cluster  at  great  distances  from  the  center ^  even  if  they  occurred  there.,  nor 
those  near  the  center  on  photographs  with  exposures  typical  of  those  used  by 
Roberts . 

We  have  examined  at  some  length  possible  correlations  for  those  12  clusters 
(which  show  obscuring  regions)  out  of  the  32  examined  by  Roberts „   Correlations 
could  be  expected  for  at  least  two  reasons s   (l)  differences  in  interstellar 
environments  and  evolutionary  activity  among  clusters 5  and  (2)  significant,, 

stematic  differences  in  the  structure  of  the  galactic  disk.  We  must  use  caution 
-  attempting  a  correlation  since  it  is  not  at  all  certain  that  the  20  clusters 
without  reported  obscuration  contain  less  diffuse  material..   Roberts  has  stated 
(private  communication.,  1963)  that  he  was  conservative  in  his  estimates  of 
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clouds  and  that  his  list  of  12  could  possibly  be  increased  by  50%. 

We  have  noted  that  only  about  half  of  the  globular  clusters  have  been  re= 
ported  as  having  RR  Lyrae  variable  stars  (Hogg,,  1959).  Adding  w  Centaur i 
(Fitzgerald  1955)  and  Kk   (idlis)  to  Roberts's  list  of  12  clusters 5  we  note 
it  RR  Lyrae  stars  have  been  sighted  in  all  but  two<>   Of  the  20  clusters 
iiined  by  Roberts  with  negative  results  for  obscuring  clouds 5  13  are  listed 
......    -laving  RR  Lyrae  stars,,   Generally  speaking 5  the  lack  of  RR  Lyrae  stars  in  a 

]  ster  is  attributed  to  the  absence  of  stars  in  the  appropriate  evolutionary 
stage.  While  the  above  figures  don't  indicate  a  signifi  ant  difference,  in 
this  case,  some  similar  criterion  having  evolutionary  significance  might,  possibly 
overed  by  a  closer  search  and  classification  of  variable  stars  and  extremely 
ue  stars  in  globular  clusters . 

1   •  sd  in  Chapter  IV  that  collisions  with  the  dense  clouds  of  the  disk 
Is  the  major  factor  for  loss  of  clouds  from  the  globulas  clusters.   The  force 
rted  by  the  inter cloud  medium  is  not  significant.  Accordingly 9   we  might 
1  that  If  there  exists  a  large s    :<   ntiguous  area  of  the  galactic  plane  with 
a  significantly  smaller  proportion  oi  dense  (10  cm  )  stn  t  ire  than  in  the 
solar  neighborhood 5  that  those  clusters  passing  through  this  region  in  their 
galactic  orbits  would  retain  more  interstellar  material  and  be  more  likely  to 
show  obscurations, 

Rougoor  and  Oort  (i960)  have  remarked  that  there  appears  to  be  a  relatively 
small  amount  of  dense  HI  structure  in  the  region  600  pc  to  3  kpc  distance  from 
the  center  in  the  galactic  plane.   Optical  observations  of  the  Sb  galaxy 

28^1  (Sandage^  1961)  show  an  amorphous  region  devoid  of  dust  or  spiral 
structure  in  the  same  region  (as  scaled  to  our  Galaxy). 

have  determined >  by  the  use  of  a  formula  from  Eggen,  Lynden-Bell,,  and 
Sandage  (1962)=  that  the  average  orbit  of  a  globular  cluster  (e  =  0.8)  with 
apogalacticon,  10  kpc.,  sweeps  out  about  205°  between  successive  apogalacticon. 
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also  have  determined  that  our  average  0,8  e  orbit  will  pass  within  3  parsecs 
bhe  center  of  the  Galaxy  during  disk  transit  if  the  semi -major  axis  lies 
within  a  cone  with  vertex  angle  of  5O0,,  centered  on  the  galactic  axis0  We 
have  attempted  a  correlation,,  without  success^,  "between  the  present  location 
of  clusters  within  the  cone  and  the  presence  of  oh s curing  material,   The  fact 
at  heliocentric  distances  of  globular  clusters  have  probable  errors  up  to  50$ 
es  any  such  correlation  doubtful,, 

Perhaps  the  two  most  significant  recommendations  for  the  future  study  of 
globular  clusters  ares   (l)  An  observational  program  that  will  determine  better 
miber  of  hot.  blue  stars  in  globular  LusterSo   For  example.,  the  incom- 
ess  factor  of  about  15  f  1  Li  M£  1  .  aid  be  improved  by  more  care- 

observations.   These  stars  are  significant  in  their  ettt    f  on   the  ionization 
d  ED  clouds a  and  the  rad  ation  field  :f  these  stars  might  be  the  most  important 
bor  influem  Lng  the  lifetime  -if  El    <  1.  ids5  si-.  -  they  have  a  reasonable  chance 
iurviving  passage  through  the  gala  T   disk.,   (2)  It  is  desirable  to  have 
a  more  detailed  program  for  observing  obscured  regions  in  globular  cluster.-;, 
primary  interest  would  be  the  photographing  ci    clusters  at  many  different 
exposures  to  determine  better  the  cloud  locations  and  their  absorption,, 
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